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ABSTRACT. 


Pulse  mode  radar  operation  Is  analyzed  under  the  assumption  that 
the  scattering  object  T  lies  in  the  far  field  of  both  the  transmitter 
and  the  receiver.  It  is  shown  that,  in  this  approximation,  the  radar 
signal  is  a  plane  wave  s(x  •  80  -  t,90)  near  T,  where  0O  is  a  unit  vector 
directed  from  the  transmitter  toward  T,  and  similarly  the  echo  is  a 
plane  wave  e(x  •  0  -  t,0,0o)  near  the  receiver,  where  0  is  a  unit  vector 
directed  from  F  toward  the  receiver.  Moreover,  it  is  shown  that 


e(T,0,0fl)  -  Re  ||  eiTu)  T(u>0,u>0o)  sfa.OoMwj 


where  s(<i),0o)  is  the  Fourier  transform  of  s(t,80)  and  T(aj0,a)0o)  is  the 
scattering  amplitude  in  the  direction  9  due  to  the  scattering  by  T  of  a 
CW  mode  plane  wave  with  frequency  u>  and  propagation  direction  9„ . 
Finally  the  singularity  expansion  method  is  used  to  show  that 


e(T,0,0o) 


iTU) 

l  e  n  Tn(6,e0)  s(u)n,0o),  Im  ion  <  0. 
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1.  INTRODUCTION  -  RADAR  ECHO  PREDICTION. 


b 


This  paper  presents  an  application  of  C.  E.  Baum's  singularity 


expansion  method  (SEX)  J2 ]  and  the  author’s  method  of  asymptotic  wave 


functions  []  1  r— L2r  13]  to  the  prediction  of  pulse  mode  radar  echoes  from 
bounded  scatterers.  The  results  presented  here  are  generalizations  of 


C 


corresponding  results  for  sonar  echoes^JrS7] .  Only  a  summary  of  the 
principal  concepts  and  results  is  presented  here.  A  complete  exposition 
of  the  theory  is  planned  for  a  separate  publication. 


1.1  Physical  Assumptions.  Radar  echo  structure  is  analyzed^J»e±trrf  under 
the  following  assumptions! 

-•  The  radar  system  (transmitter  and  receiver)  operates  in  a 
stationary  homogeneous  isotropic  unlimited  medium, 

The  system  is  stationary  with  respect  to  the  medium, 

•  The  scatterers  are  bounded  perfectly  conducting  objects, 

•  The  scatterers  are  stationary  with  respect  to  the  medium. 

•  The  transmitter  and  receiver  are  in  the  far  field  of  the 
scattering  objects. 

In  addition  it  is  assumed  that  secondary  echoes  due  to  the  radar  system 
components  are  negligible.  -/ 


1.2  Mathematical  Formulation.  A  fixed  Cartesian  coordinate  system  is 
used  throughout  the  paper,  x  -  (xltx2,x3)  €  R3  denotes  a  coordinate 
triple  of  this  system  and  t  €  R  denotes  a  time  coordinate.  T  denotes  a 
closed  bounded  subset  of  R3  that  represents  the  scatterers  and  H  *  R3  -  T 
denotes  the  domain  exterior  to  T.  The  common  frontier  of  T  and  0,  which 


represents  the  surface  of  the  scatterers,  is  denoted  by  aft.  The  medium 
filling  ft  is  characterized  by  a  dielectric  constant  e  and  a  magnetic 
permeability  y.  It  will  be  assumed  that  e  =  1  and  y  -  1  since  this  can 
be  achieved  by  a  suitable  choice  of  units. 

The  electric  and  magnetic  fields  will  be  described  by  their 
components,  (E^Ej.Ej)  and  (H1,H2,HJ)  respectively,  relative  to  the 
fixed  Cartesian  system.  It  will  be  convenient  to  use  the  notation  and 
conventions  of  matrix  algebra  and  to  characterize  the  electromagnetic 
field  by  the  6*1  column  matrix 


(1.1) 


u  -  u(t,x)  -  (Ej  E2  E3  Hj  H2  H3) 


where  M  denotes  the  transpose  of  matrix  M.  Similarly,  if  the  electric 
and  magnetic  current  densities  that  generate  the  field  are  described  by 
their  components,  (Jj.Jj.Jj)  and  (Jj,J’,JJ)  respectively,  then 


(1.2) 


f  -  f(t,x)  -  (Jj  j2  j3  Jj  jpT 


characterizes  the  field  sources.  With  these  conventions  Maxwell's 
field  equations  can  be  written 


(1.3) 


)  u  +  £  A,  D.u  +  f  -  0  for  t  S  R,  x  €  ft 

C  j-1  J  J 


where  Dt  -  9/3t,  Dj  »  9/3xj  (j  ■  1,2,3)  and  At,  A2,  A3  are  the  three 
symmetric  6  *  6  matrices  defined  by 


3  f  0  «(p,l  °  P>  'P'' 

(1.4)  l  A,  p.  -  ,  M(p)  -  -p  0  p 

J-1  J  J  l-M(p)  0  1 

P2  "Pi  0  . 
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The  field  equations  (1.3)  will  be  supplemented  by  the  boundary  condition 
for  a  perfect  electrical  conductor.  It  can  be  written 

(1.5)  M(n)  E  -  0  on  3Q 

T 

where  n  **  (nltn2,n3)  is  a  unit  vector  on  3Q  and  E  «  (Ej  E2  E3)  is  the 
electric  part  of  u. 

A  theory  of  solutions  with  finite  energy  of  (1.3),  (1.5)  was 
given  in  [10] .  The  total  field  energy  at  time  t  is  given  by 


(1.6) 


E  “  J  J  u(t,x)T  u(t,x)  dx 


where  dx  «  dx,  dx2  dx3 .  The  theory  of  [10]  makes  use  of  the  energy  norm 


(1.7) 


lulfl  *  U(X)T  u(x)  dxj 


,1/2 


and  corresponding  Hilbert  space  3f.  The  pulse  mode  radar  echoes 
constructed  below  are  in  3f. 


> 


» 


♦ 


» 


I 
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2.  PULSE  MODE  RADAR  SIGNAL  STRUCTURE. 

The  transmitter  will  be  assumed  to  be  localized  in  the  ball 
B(x0,6fl)  »  (x  :  |x  -  x0J  <  S0)  and  to  act  during  an  interval  0  <  t  <  t0. 
Thus  the  source  distribution  f  in  (1.3)  will  have  support 

(2.1)  supp  f  C  {(t,x)  :  0  <  t  f  t0  and  |x  -  x#  |  <  <S0}  . 

The  corresponding  pulse  mode  radar  signal  is  the  electromagnetic  field 
u0(t,x)  that  is  generated  by  f  when  no  scatterers  are  present.  Thus  uQ 
is  characterized  by  the  conditions 

3 

(2.2)  D  uj  +  l  A,  D ,u0  +  f  -  0  for  t  £  R,  x  6  R3, 

j-1  J  J 

(2.3)  u0(t,x)  «  0  for  t  <  0,  x  €  R3. 

The  field  u0  can  be  constructed  by  Fourier  analysis  or  by  the  method  of 
retarded  potentials  [11,  13,  15]  but  these  constructions  will  not  be 
recorded  here. 

2.1  Asymptotic  Wave  Fields.  For  definiteness  the  scatterers  are 
assumed  to  be  localized  in  the  ball  B(0,5)  centered  on  the  origin: 

T  C  B(0,6).  With  this  convention  the  assumption  that  the  transmitter 
lies  in  the  far  field  of  T  can  be  formulated  as  |x0|  »  1.  The  signal, 
propagating  at  the  speed  c  ■  (eu)”1^2  ■  1,  will  arrive  at  F  at  a  time  t 
of  the  same  magnitude  as  [xg|,  whence  t  >>  1. 

It  was  shown  in  [13]  that  each  signal  u„  with  finite  energy  has 
an  asymptotic  wave  field  u*  of  the  form 


J 
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(2-4)  u0(t,x)  -  |x  -  x0|  1  s(|x  -  x0 |  -  t,0),  0  =  (x  -  x0)/|x  -  xQ | 

such  that 

(2.5)  lim  I u„ (t, •)  -  u“(t,*)lR3  =  0  . 

t**+°° 


The  wave  profile  s(x,9)  is  defined  for  all  (T,0)  £  R  x  s2  where  S2  is 
the  unit  sphere  in  R3.  Moreover,  by  specializing  the  results  of  [13]  it 
can  be  shown  that  s(x,0)  has  the  properties 


(2.6) 


s(x,0)T  s(x,0)  d©  dx  <  00  , 

R  >  S2 


where  d0  is  the  element  of  area  on  S2  (solid  angle) ,  and 
(2.7)  P(0)  s(x,0)  -  s(x,0) 


where 


(2.8) 


P(0) 


1  -  00 
-M(0) 


M(0) 

1  -  00, 


for  all  0  6  S2. 


In  (2.8),  00  denotes  the  dyadic,  or  tensor,  product  of  0  with  itself 
with  components  ©j 9^ .  Property  (2.7),  (2.8)  characterizes  the  polariza- 

CO 

tion  properties  of  the  asymptotic  wave  fields  uQ . 

The  function  s(x,9)  will  be  called  the  pulse  mode  transmitter 
radiation  pattern.  It  can  be  constructed  from  the  source  function  f; 
see  [15].  However,  it  will  be  assumed  here  that  s,  rather  than  f,  is 
given  since  s  is  the  important  function  in  pulse  mode  transmitter 
design.  The  construction  of  a  transmitter  with  a  prescribed  radiation 
pattern  is  the  task  of  the  transmitter  design  engineer. 
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2.2  The  Plane  Wave  Approximation.  Define  0O  e  S2  by  xQ  -  — | x0  J  90 . 
Then  0O  is  directed  from  the  transmitter  toward  the  scatterers  and  for 
x  near  r  one  has 

(2.9)  |x  -  x0|  =  |  x0 1  +  0O  •  x  +  0(|xj_1)  for  |x0|  »  1. 

Hence,  by  (2.4), 

(2.10)  u>) » i*0ri  s(0o  •  x  - 1  +  |xoi,0o)  +  o(|x0r2) 

near  F.  If  the  error  term  is  dropped  one  has  a  pulse  mode  plane  wave 
signal.  This  approximation  is  made  in  the  remainder  of  the  paper. 


3.  PULSE  MODE  PLANE  WAVE  SCATTERING. 

A  plane  wave  signal 

(3.1)  u0(t,x)  »  s(x  •  60  -  t, 90) ,  supp  s(*,0Q)  C  [a,b] , 
is  assumed  where  the  wave  profile  s(t,9o)  satisfies 

(3.2)  P(0O)  a(T,e0)  =  s(t, 90)  . 

Such  a  field  is  a  solution  of  Maxwell's  equations  (2.2)  with  f  =  0.  The 
total  field  u(t,x)  resulting  from  the  interaction  of  u0(t,x)  with  the 
scatterers  is  characterized  by  the  properties 

(3.3)  Du  +  £  A.  D.u  ■  0  for  t  6  R,  x  e 

j-1  J  3 

(3.4)  M(n)  E  =  0  for  t  e  R,  x  6  3(2, 

(3.5)  u(t,x)  *  u0(t,x)  for  t+b+S<0,  x£f! 

where  E  *  (ut  u2  u3)  is  the  electric  part  of  u.  The  scattered  field, 
or  echo,  is  defined  by 

(3.6)  ug(t,x)  *»  u(t,x)  -  u0(t,x)  for  t  e  R,  x  6  Q. 

The  author  has  shown,  by  the  method  of  [11,  15],  that  ug  has  an 
asymptotic  wave  field 

(3.7)  u“(t,x)  -  |xj-1  e ( 1 x j  -  t,e,0o),  x  -  jxj  9 


9 
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that  converges  to  ug(t,x)  in  energy  when  t  -*•  °°: 

(3.8)  lim  Sue(t,*)  -  u“(t,*)Uo  =  °* 

t"*+#0 

The  proof  follows  that  for  the  scalar  case  of  [15]. 

Points  x  in  the  far  field  of  T  satisfy  jx|  »  1.  The  echo  u& 
will  arrive  at  a  receiver  at  such  a  point  when  t  »  1.  Hence  the  echo 
may  be  approximated  in  the  far  field  by  the  asymptotic  field  (3.7).  For 
this  reason  e(T,0,0o)  will  be  called  the  echo  waveform.  It  depends  on 
the  direction  of  incidence  of  the  plane  wave  (3.1)  and  the  direction  of 
observation  0.  In  this  approximation,  the  echo  prediction  problem  is 
the  problem  of  constructing  e(t,0,0Q)  when  the  transmitter  radiation 
pattern  s(T,0o)  and  the  scatterers  T  are  given.  The  solution  to  this 
problem  given  below  is  based  on  the  theory  of  CW  mode  radar  echoes 


outlined  in  the  next  two  sections. 


4.  CW  MODE  SIGNAL  STRUCTURE. 

The  CW  mode  electromagnetic  fields  are  solutions  of  the  field 
equations  (1.3)  of  the  form 

(4.1)  u(t,x)  =  e-1Wt  v(x),  f(t,x)  =  e~i0Jt  p(x) 

whence 


3 

(4.2)  £  A.  D.v  -  icjv  =  p. 

J-l  3  3 


CW  mode  signals  i”  R 


(4.3)  G(x,x',u)  ■ 


are  generated  by  the  Green's  matrix 


W  +  a)2  13  -iutt(V) 


ioJM(V)  VV  +  w2  13 


eiw | x-x ' | 
4tto)  I  x-x '  I 


[7] 


where  1  denotes  the  n  x  n  unit  matrix.  G  is  the  outgoing  solution  of 
n 

the  equation 


(4.4) 


f  I  A.  D.  -  iu)|  G(x,x',u>)  =  6 (x  -  x')  1  . 
tj-1  J  3  J 


The  outgoing  solution  in  R3  of  (4.2)  is 


(4.5) 


v(x) 


G(x,x,,cj)  p(x')  dx' 
R3 


Asymptotic  evaluation  of  v(x)  for  large  |x|  using  (4.3)  and  (4.5)  gives 


the  far  field  form 
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(4.6)  v (x)  =  (2it)  1/2  uj  — -  P (6)  p(-o)6)  +  0(|x|~2) 

1*1 

where  x  =  |x|9,  P(9)  is  defined  by  (2.8)  and 

(4*7)  P(p)  =  ^2tt)  3/z  {  3  e_iP  X  P(x)  dx 

is  the  Fourier  transform  of  p(x) .  In  particular,  noting  that 
P(-9)  P(9)  =  0,  it  is  seen  that  the  Silver-Miiller  radiation  condition 
for  v(x)  can  be  written 

(4.8)  P(-9)  v(  |  x|  0)  =  0(|x|-2),  |  x  |  +». 

4.1  CW  Mode  Plane  Waves.  G(x,x',u))  represents  a  CW  spherical  wave  from 
a  point  source  at  the  point  x'.  On  putting  x'  =  -|x'|n  in  (4.3)  and 
making  |x'|  -+  °°  with  x  fixed  one  finds  after  a  short  calculation 

(4.9)  G(x,x',uO  =  ( 2-rr | x ’  1 ) “ 1  «  eia)lx'l  eiu)n‘x  P(n)  +  Odx'l"2). 

Dropping  the  error  term  gives  a  matrix  CW  mode  plane  wave  electromagnetic 
field.  The  general  CW  mode  plane  wave  field  is  obtained  by  applying 

(4.9)  to  a  constant  vector  and  dropping  the  error  term.  It  has  the  form 

(4.10)  v(x)  =  eip  x  P(n)c,  p  =  | p ( f| 

where  c  is  an  arbitrary  6-component  vector.  This  may  also  be  derived 
from  (3.1),  (3.2)  by  taking  s(T,n)  •  eiwT  P(ri)c.  (4.10)  is  equivalent 
to  the  familiar  formulas 


5.  CW  MODE  ECHO  STRUCTURE. 


The  columns  of  the  6  *  6  matrix-valued  function 

(5.1)  ^(x.p)  -  (2it)“3/2  eip*x  p(n),  p  -  Ip  In, 

are  CW  mode  plane  waves  of  the  form  (4.10).  The  scattering  of  the  CW 
mode  matrix  plane  wave  (5.1)  by  T  produces  a  CW  mode  matrix-valued  field 

(5.2)  'f(x.p)  -  ^(x.p)  +  ^(x.p),  x  €  ft,  p  €  R3  -  (0} 


that  is  characterized  by  the  properties 


(5.3) 


I  A.  D 


J-l 


J  3 


* (x.p) 


0,  x  6  £2, 


(5.4)  M(n)  fgU.p)  «  o,  x  e  3ft 

(5.5)  P(-0)  't'SC(|x|0,p)  -  0(|x|"2),  | x |  -*■  », 


where  is  the  electric  part  of  f  (a  3  x  6  matrix) .  The  author  has 

Ci 

shown  the  existence  and  uniqueness  of  ¥(x,p)  for  a  large  class  of  domains 
ft,  including  the  "cone  domains"  of  N.  Week  [9]  and  domains  having 
S.  Agmon's  "restricted  cone  property"  [1].  The  proofs,  which  generalize 
the  results  of  [11]  to  Maxwell's  equations,  are  based  on  compactness 
results  of  N.  Week  [9]  and  C.  Weber  [8],  respectively.  In  the  special 
case  that  3ft  is  a  smooth  surface  ¥(x,p)  can  be  constructed  by  the 
integral  equation  method  described  below. 

sa 

5.1  Far  Field  Form  of  CW  Mode  Echoes.  ¥  (x, p)  is  the  CW  mode  echo 

produced  by  the  scattering  of  'i|0(x,p)  by  T.  An  integral  representation 
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of  'i'sc  by  the  Green's  matrix  (4.3)  can  be  used  to  derive  the  far  field 
form 

i  |  p | |x| 

(5.6)  fsc(x,P)  =  ^ - Klpje.p)  +  C7(|x|-2).  x  =  |x|e, 

4tt  |x  | 

where  T(p,p')  is  a  6  x  6  matrix-valued  scattering  amplitude.  The 
polarization  of  the  echo  in  the  far  field  is  characterized  by  the 
property 

(5.7)  P(n)  T(|P|ri,  | p|n ')  =  0. 


5.2  Construction  of  T(p,p').  Define 

(5.8)  J(x,p)  *  n(x)  x  ^(x,?),  x  6  9(2, 

where  ^  is  the  magnetic  part  of  ¥.  J(x,p)  is  the  matrix  electric  current 
density  on  9(2  induced  by  the  plane  wave  0 .  The  divergence  theorem  and 
the  jump  relations  of  potential  theory  can  be  used  to  show  that 


(5 


.9)  J(x,p)  »  2(n  x  ¥„(x,p))  +  [  K(x,x',|p|)  J(x',p)  dS 

'9(2 


where  K  is  the  3x3  matrix-valued  kernel 


(5.10)  K(x,x',uj) 


2tt 


iui|x-x’j 

V  - -  n(x) 

|x-x'| 


l.e1"l«,l  ,  1 

9n 


|x-x' 


If  9(2  is  smooth  then  (5.9)  is  a  Fredholm  equation  and  can  be  used  to 
construct  J(x,p)  and  'f(x.p);  cf.  L.  Marin  and  R.  W.  Latham  [4]  and 
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L.  Marin  [5].  The  scattering  amplitude  can  be  calculated  from  J(x,p) 
and  the  relation 


(5. 


1  ^“(x.p)* 

J(x.p’) 

ha  1 

.  0  , 

ds,  |p|  -  ]p'| 


! 


6.  PULSE  MODE  RADAR  ECHO  STRUCTURE. 

The  solution  of  the  pulse  node  radar  echo  prediction  problem 
formulated  in  §3  is  given  by  the  relation 

(6.1)  e(T,6,80)  -  Re 

where 

(6.2)  sCMo)  -  -(Yffprr-  I"  e~iU)T  s(T,e0)  dT 

is  the  Fourier  transform  of  s(T,60).  Thus  under  the  far  field 
assumptions  of  Si  the  echo  waveform  is  determined  by  the  transmitter 
waveform  and  the  matrix  scattering  amplitude  T(ti>0,a>0o) .  The  latter  can 
be  calculated  by  solving  the  integral  equation  (5.9)  and  using 
relation  (5.11). 

Equation  (6.1)  is  the  generalization  to  electromagnetic  fields 
of  the  analogous  result  for  acoustic  scattering  that  was  derived  in  [15]. 
A  proof  of  (6.1)  may  be  given  by  the  method  of  [15].  The  key  item  in 
the  proof  is  the  theorem  that  the  CW  mode  fields  ^(x.p)  are  a  complete 
family  of  generalized  eigenfunctions  for  the  Maxwell  system.  A  proof 
along  the  lines  of  [11]  may  be  based  on  the  results  of  Week  [9]  or 
Weber  [8]. 


(f 


e*TU  T(a>8,u)0o)  s(u>,9o)  du 
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7.  SEM  EXPANSION  OF  PULSE  MODE  RADAR  ECHOES. 

If  Che  scatterers  T  are  bounded  by  smooch  surfaces  Che  incegral 
equation  (5.9)  can  be  solved  for  J(x,io0)  by  Che  Fredholm  decerminanc 
mechod  [14].  Note  thac  to(x,<o0)  and  K(x,x',<o)  are  entire  functions  of 
<o.  It  follows  from  Che  Fredholm  theory  Chat 


(7.1) 


and  hence 


D(co) 


(7.2) 


T(u>0,io0..)  ■ 


where  D(<o),  M(x,<o8)  and  N((o6,<o0o)  are  entire  functions  of  to.  Moreover, 
the  poles  of  T((O0,<o0o)  can  be  shown  to  lie  in  the  lower  half-plane. 
These  facts  can  be  used  to  develop  an  SEM  expansion  of  the  echo 
waveform  (6.1). 

The  reality  of  s(x,90)  and  symmetry  properties  of  T(p,p')  imply 
that  (6.1)  can  be  rewritten 


(7.3) 


:,0.0«)  -  T  f 

*  -oc 


eiTU  T(!o0,io0o)  S(u,e,)du. 


It  is  natural  to  regard  this  integral  as  a  contour  integral  in  the 
arplane  and  to  shift  the  contour  to  a  line  Im  to  ■  -b  <  0.  Assume  that 
the  poles  <oQ  of  T(<o0,io0o)  satisfy 


(7.4) 


D'(ion)  +  0,  n  -  1,  2,  3,»** 
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(7.5) 


{n  :  -b  <  Im  oj  <  0}  is  finite 
■“  n 


(7.6)  |N(u)0,weo)  |  <  C|ui|m  for  -b  <  Im  u  <  0 

vhere  C  and  m  are  constants.  Then  (7.3)  implies 


(7.8)  Tn(6,60)  -  -iri  Res  T(aj6,u)90). 

1 

Hypothesis  (7.4)  is  inessential.  If  T Coo0 , to0 0 )  has  a  higher  order  pole 

it^n 

then  in  (7.7)  e  will  be  multiplied  by  a  polynomial  in  T.  Hypotheses 
(7.5)  and  (7.6)  are  closely  connected  with  the  geometry  of  T  and  the 
associated  question  of  the  exponential  decay  on  bounded  sets  of  the 
scattered  fields.  For  acoustic  scattering  there  is  a  considerable 
literature  on  these  questions:  see  [3]  and  [6]  and  the  literature  cited 
there.  The  electromagnetic  case  awaits  further  analysis. 
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ABSTRACT. 

Pulse  mod*  radar  operation  is  analyzed  under  Che  assumption  chat 
eh*  scattering  object  T  lies  in  the  far  field  of  both  the  transmitter 
and  the  receiver.  It  is  shown  that,  in  this  approximation,  the  radar 
signal  is  a  plane  wave  six-  9Q  -  c,90)  near  T,  where  6fl  is  a  unit  vector 
directed  from  the  transmitter  coward  ",  and  similarly  the  echo  is  a 
plane  wave  e(x  •  9  -  t,9,0j)  near  the  receiver,  where  9  is  a  unit  vector 
directed  from  f  toward  che  receiver.  Moreover,  it  is  shown  that 

«<r,e.sj)  -  R«  ||  •  tt“  Tt-e.-e,)  »(w,5,)du,| 


where  8(01,9,)  is  the  Fourier  transform  of  s(t.9,)  and  T(wd,u.d3>  is  che 
scattering  amplitude  in  che  direction  a  due  to  che  scattering  by  f  of  a 
CM  mode  plane  wav*  with  frequency  w  and  propagation  direction  9a. 
Finally  the  singularity  expansion  method  Is  used  to  show  that 


(T.3,9,)  '  l 
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